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ABSTRACT 


In April 1969, 50 different gravity stations on and around the 
perimeter of Fort Ord, California, were obtained using a LaCoste- 
Romberg Model G Geodetic Gravity Meter. The density of stations en- 
abled accurate location of 5-mgal contours of Simple Bouguer Anomaly. 
The major differences found between the new and previously published 
contours occurred in the west and southwest regions of Fort Ord. 

Five stations obtained by an earlier investigator were re- 
occupied during this study. The differences in observed gravity at 
these stations ranged from -0.37 mgal to -0.70 mgal. 

In an attempt to explain the differences, 1] additional stations 
were reoccupied in August 1969. Observed gravity differences for 
these stations ranged from -0.05 mgal to -0.58 mgal. The differences 
could not be fully explained nor could the earlier study be success- 
fully tied to this study. 

Fort Ord lies on a gravity low and is isostatically overcompen- 
Sated. Further gravity readings are required on the Monterey 
Peninsula and in the Salinas Valley to adequately define the sub- 


structure of Fort Ord. 
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I INTRODUCTION 

Gravity measurements are made for a number of reasons. Some of these 
are scientific and some are commercial. The geophysicist uses gravity 
anomalies as an aid in interpreting the geological substructure of an 
area. The mining or petroleum geologist 1s looking particularly for 
promising structures for commercial exploitation. Geodesists use 
gravity anomalies as an aid in determining the shape of the earth. 
Gravity values, once ascertained, may be of future use to the Armed 
Forces as a means of precise positioning. 

The original purposes of this study were (1) to obtain practical 
experience in conducting a complete gravity survey of an area, (2) 
to reduce the data obtained to a standard usable form, and (3) to de- 
termine the geological substructure of the area surveyed. 

The first two objectives were successfully completed. The third 
purpose was not successfully concluded for the reasons explained in 


the final pages of this paper. 


Dy, re SWRVEY 


A. SELECTION OF AREA 

An unsurveyed area close to Monterey, California, was desired so 
that the survey could be conducted during weekends and half-days that 
were available. Two possible areas met the above criteria: (1) the 
northern end of the Santa Lucia Range located to the south of the 
Monterey Peninsula and (2) Fort Ord, California. The former area is 
roadless in the unsurveyed areas. Neither the time nor the resources 


required to conduct a horseback survey of the area was available to 


the author. The latter and smaller of the two areas was therefore chosen. 
Fort Ord has roads covering most of its previously unsurveyed area and 
its proximity to the Naval Postgraduate School made three-hour field 

EV lps migaet lca ll. 

Previous gravity surveys around the perimeter of Fort Ord indicated 
the existence of a local low Simple Bouguer Anomaly [ Bishop and Chapman, 
1967] , but no observations had been made on Fort Ord to accurately lo- 
cate the site of the low. Arrangements were made by the author to 
enter the firing range area to obtain stations which would help deline- 


ate this low. 


B. EQUIPMENT 

Two gravimeters were made available to the author by the U.S. 
Geological Survey (U.S.G.S.). A LaCoste and Romberg Model G. Geodetic 
Gravity Meter, serial G58, was utilized for all observations made in 
April 1969. Serial number G143 was used for the August 1969 observa- 
tions. 

Both meters have a range greater than 7000 milligals (mgal), a 
reading accuracy of £ 0.01 mgal, and a drift rate of less than 1.0 mgal 
per month. Model G meters are sealed to eliminate any effect from 
changes in atmospheric pressure and, as a safety precaution, are also 
internally pressure compensated. The gravity sensor is completely de- 
magnetized and enclosed with a magnetic shield. The entire meter is 
maintained at a constant temperature while in operation. 

A simplified diagram of the basic meter is shown in Fig. |. The 
gravity response system consists of a weight on the end of a horizon- 
tal beam supported by a zero-length spring. The shock-eliminating 


springs form a floating pivot, thus eliminating any friction in the 


Connecting Links 


lever 
: Measuring Screw 






Shock Eliminating Spring 


Fig. |. Simplified Diagram of Meter 
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moving system. The gravity response system is completely suspended by 
Springs and will withstand any shock that does not damage its housing. 

The lever system and measuring screw shown in Fig. 1 are accurate 
ly calibrated over their entire range. Since the calibration factors 
depend only on the quality of the measuring screw and the lever system 
and not upon any type of auxiliary springs they do not change per- 
ceptibly with time [LaCoste- Romberg | ‘ 

Topographic maps at a scale of 1:24000 of the Fort Ord area were 
provided by U.S.G.S. to provide elevation control for the observations. 
Additional data was provided by the U.S. Army Corps of Engineers from 
their elevation control surveys conducted by the 663 Engineering Com- 
pany in April 1962 and by the Sacramento District Corps of Engineers 


in February and March 1969. 


C. DATA COLLECTION 

During the period 3-13 April 1969, 51 different stations were 
occupied on six separate days. For each set of runs the first and 
last station occupied was the California Division of Mines and Geo- 
logy Base Station No. 259 (CH259) located at the Monterey Airport. 

Of the 51 stations, five were obtained and reported earlier 
[Sieck, 1964) ,» one was the base station and the remaining 44 
stations were new. 

An attempt was made to locate and occupy all monumented bench 
marks around the perimeter of or on Fort Ord which were shown on 
the U.S.G.S. topographic maps. Six of 15 were recovered and oc- 
cupied. The remaining nine monumented bench marks were not found 
by the author. A few of these had obviously been destroyed by con- 


struction and others were most likely concealed by thick vegetation. 


The remaining stations were located at non-monumented bench marks 
shown on the U.S.G.S. topographic maps, monumented bench marks, both 
permanent and temporary, established by the U.S. Army Corps of Eng- 
ineers and, in one case at the peak of a hill. This last station was 
the only one where the elevation was based on a spot elevation rather 
than a vertical control survey. 

At each station occupied the station number, meter reading, time 
of reading, height relative to the reference elevation, reference 
elevation, and description of the surrounding terrain were recorded. 

At the conclusion of each set of runs the two readings taken at 
the base station were corrected for tide effects and compared in 
order to check for tares and to compute drift corrections. The 
drift was assumed to be linear over the period of the observations. 
The maximum drift rate used was found to be 0.04 mgal/hr. This drift 
rate was applied over two stations on 3 August. During the original 
Survey in April the drift rates ranged from 0.00 mgal/hr to 0.02 
mgal/hr. 

One set of observations on 3 August was discarded due to an ap- 
parent tare of about 0.25 mgal, since the drift rate required to 
account for the change in observed values at the base station would 


have been 0.10 mgal/hr. 


ITT. GRAVITY CORRECTIONS 


Gravity observations must be corrected to reduce the measurements 


to the value which would be observed on a uniform spheroid fitted as 


Closely as possible to sea level. This reference spheroid is competely 
smooth with the vertical distribution of density with depth horizontally 


uniform. 


A. THEORETICAL GRAVITY 

The reference spheroid is best approximated by a triaxial spheroid, 
which is an ellipsoid of revolution modified by depressions along the 
tvo 45°_latitude lines and by a flattening and bulging of the equator 
[ Dobrin, 1960 | . The generai form of the gravity variation to be 
expected on such a reference spheroid is given by: 


ane Cans Cp sient (1) 


where g, is the value of gravity at the equator at longitude 180°, 


C, and C, are constants which give a measure of the earth's true shape 


/ 
and L is latitude: 

The constants are determined by fitting world-wide pendulum measure- 
ments of gravity, which have been reduced to sea level, to the formula 
by the method of least squares. Since 1884 various investigators have 
determined the value of the constants using measurements available 
at the time. 

The 1930 International constants were used for this study, which 
when substituted into (1) gives the theoretical gravity (THG) in 
mgal at any latitude L: 


THG = 978049 [ 1 + 0.0052884 sin* L - 0.0000059 sin* 2L] (2) 


B. TIDAL CORRECTION 
Before the actual gravity values as read at each station were re- 
duced to sea level the effect of the orbiting moon was removed. The 


orbital parameters of the moon are well known and the tidal effect 
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can be predicted in advance. The tidal gravity corrections were de- 
termined by use of the tables and nomogram of the Service Hydrographique 
de la Marine and Compagnie Generale de Geophysique [1968 | : 

The tidal correction and instrument drift correction were combined 
and added to the measured value of gravity to obtain the observed 


gravity. 


C. FREE AIR CORRECTION 

The next correction that was made to the measured gravity was the 
free air correction (FAC). This correction assumes the existence of 
Only air between the station and sea level and accounts only for the 
difference in elevation. 

The correction in mgal applied for this survey was: 

FAC = [ 0.09411549 - 0.000137789 sin* L J E 
0. 0000000067E (3) 

where E = elevation in feet of station above sea level. 
| Scheibe and Howard, 1964 | (See Appendix 1 for a derivation of this 


correction. ) 


D. BOUGUER CORRECTION 
An adjustment must be made for the oversimplification of assuming 

Only air exists between the station and sea level. This was done by 
assuming that a slab of density P > thickness E and extending infinite- 
ly in all horizontal directions exists between the station and sea 
level. The effect of this slab is called the Bouguer Gm@mrection 
(BC) and is given in mgal by: 

BC = 27/¥ gE= 0.012774 c E (4) 
where fs ecens icy sin g/cm? 

¥ = univeral gravitational constant 


[ Grant and West, 1965 and S. L. Robbins, personal communication | 
ie 


FE. TERRAIN CORRECTION 

Another correction was required to account for the fact that there 
is not a uniform slab, whose upper surface is the "Bouguer plane", 
but rather there is the actual topography of the land in the vicinity 
of the station. Both the material extending above the hypothetical 
Slab and the material missing from the slab due to the actual topo- 
graphy result in a reduction of the measured gravity values. Hills 
reduce the value because the mass opposes the attraction at the 
Station due to the mass of the earth, and the valleys also cause a 
reduction because they were included in the Bouguer correction, but 
in fact are not there. These corrections are added to the observed 
gravity. 

In theory the volume and average density of each variation from 
the Bouguer plane must be calculated and its effect on the station as 
a function of distance must be determined. In practice a template is 
used which divides the area in the vicinity of the station into com- 
partments. The template consists of concentric circles with each 
annular ring divided into compartments. In general the number of com- 
partments in each ring increases with increasing radius. The average 
elevation in each compartment is determined by visual inspection. 

The absolute value of the difference in elevation between the com- 
partment and the station is used to enter a correction table which 
provides the applicable correction for a given assumed density. 

Hayford-Bowie templates were used throughout this survey through 
Zone F. For terrain corrections at distances between 2.29 km and 
166.70 km from the station, the U.S.G.S. computerized Terrain Cor- 


rection Program was utilized. 


ie 


The principle features of this program are that topography is digi- 
tized on a latitude and longitude grid, compartments are assembled on 
the basis of map boundaries and an exact tie to a circular inner 
boundary (the outer radius of zone F) is made. 

All terrain corrections (TC) were based on an assumed density of 


2.67 gm/cm- 


F. CURVATURE CORRECTION 

A final correction was required to account for the curvature of the 
earth. The Bouguer correction assumes a flat earth which is a reason- 
able assumption for short distances, but is inaccurate for the greater 
distances utilized by the computerized Terrain Correction Program. 

Lambert [1930] developed the curvature correction (CC) to account 
for the fact that the earth is not flat. Although he computed the 
correction to four significant places in the first term, he utilized 
the value of 77 only to three significant niece U.S.G.S. has recom- 
puted the correction [ S. L. Robbins, personal communication ] and 
found it to be in mgal: 


-8§ 92 —/S 3 
CC = 0°0004462E -- 3.252 x WOP bk Sic x 10 e (5) 


IV. GRAVITY ANOMALIES 


A gravity anomaly is the residue left after corrections have been 
applied to the observed gravity and the theoretical gravity has been 
subtracted from the result. There are different names for gravity 


anomalies, depending upon which corrections have been applied. 


A. FREE AIR ANOMALY (FAA) 
This anomaly is given by: 
FAA = 0G + FAC - THG | (6) 
B. SIMPLE BOUGUER ANOMALY (SBA) 
The SBA is given by: 
SBAe= 0G +.FAC =9BC-oaiG (70) 
It is this anomaly that is normally used for comparative purposes 


and. thateswhiehpis plotted on Fig. 3 aombigaed/. 


C. COMPLETE BOUGUER ANOMALY (CBA) 
When all corrections discussed above are included the result is 
the Complete Bouguer Anomaly: 


CBA = 0G + FAC - BC + TC - CC - THG (8) 


V. DATA REDUCTION 


Manual corrections for tide, drift and terrain corrections through 
Hayford-Bowie Zone F were computed by the author. All other data re- 
duction was accomplished by existing U.S.G.S. Gravity Reduction Pro- 
grams on the U.S.G.S. IBM 360 Computer. The results are summarized 


beginning on page 48 . 


ViteeRESULTS 


As shown in Fig. 2 to 27 a local Simple Bouguer Gravity low of 
-35 mgal exists in the extreme eastern portion of Fort Ord. This 
low decreases to the southwest, becoming zero at approximately the 


Monterey - Seaside line and continuing positive to the southwest. 
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On Fort Ord itself the SBA ranges from - 5.05 mgal to -35.59 mgal. 
A value of - 38.52 mgal was located on Reservation Road which runs 
between Salinas and Marina forming the eastern boundary of Fort Ord. 

The large number of stations in the relatively small area allowed 
the 5-mgal contours of SBA to be accurately located. The new contours 
agree well with the contours plotted by Bishop and Chapman [1967 | 
in the eastern portion of Fort Ord. 

In the western and southwestern area there are major changes in 
the contour location. These changes are attributed to the greater 
saturation of stations made during this study. 

During April 1969, five stations were recovered which had been 
obtained previously and reported by Sieck [1964]. Table I compares 
the values of observed gravity at these five stations. 

An attempt was made to account for the observed differences, but 
no satisfactory explanation was evident. The first and most obvious 
explanation was that CA 879, CA 883 and CA 843 were obtained on one 
day and the remaining two stations were measured on another day. Ivey's 
data shows that CA 890 and CA 843 were surveyed on April 6, 1969, and 
the remaining three stations were occupied on April 3, 1969. Sieck's 
data does not include dates, but the numbers he assigned to the 
original stations suggest that CA 884 and CA 883 were surveyed on the 
Same day, CA 879 and CA 890 on some other day and CA 843 on still 
another day. No correlation between observed gravity and the date 
of the measurement could be found. 

A portion of the observed differences may possibly be attributed 
to the fact that Sieck used Worden gravimeters provided by U.S.G.S. 


These meters can not be read to better than + 0.1 mgal and a figure 


of + 0.15 mgal is probably more realistic. 


four base stations in the area for his use, and tied them to CH 176 


Further, Seick established 


(Woollard Station WU-3) located at Stanford University. 


COMPARISON OF OBSERVED GRAVITY MEASURED AT FIVE STATIONS 


OBSERVED GRAVITY 
(+979800 mgal ) 


Station STECK ay Difference 

CA 879 62.92 63.40 -0.48 

CA 884 40.07 40.81 -0.64 

CA 883 49.73 Oe -0.38 

CA 890 43.05 43.75 -0.70 

CA 843 33.50 SOmOT =0no7 
TABLE 1 


The calibration error for the U.S.G.S. Worden gravimeters has 
been determined to be one in six hundred. Thus for every six mgal 
difference in gravity between two stations the meter will introduce 
an error of 0.01 mgal. The difference in gravity between Stanford 
and Fort Ord is about 78 mgal which would produce an error of 
approximately 0.13 mgal. 

If both of the above errors are additive a possible error of 
about 0.33 mgal may exist in Sieck's data. To confirm this, during 
the period 1 - 3 August 1969, an additional 1] stations occupied 
by Sieck were recovered by Ivey. These 11 stations were those in 
Marina, Seaside, Spreckles and Salinas (west of US 101) Quad- 


rangles that could be found. Table II shows a comparison 


of observed gravity as measured at these I] stations. If the same errors 
are assumed to apply to Seick's data the differences range from + 0.38 
mgal to - 0.25 mgal. 
COMPARISON OF OBSERVED GRAVITY AT 
1] DIFFERENT STATIONS IN MONTEREY-SALINAS VICINITY 


OBSERVED GRAVITY 
( + 979800 mgal) 


Station STECK IVEY Difference 
CA 845 Jen OU Geo O05 
CA 844 62.20 62.27 -0.07 
CA 716 D1 oO) S705 -0.13 
CA 873 58. 40 98.56 -0.16 
CA 702 3070 99.88 -0.18 
CA 886 Zon © Zon 90 -0.22 
CA 709 45.70 46.06 -0. 36 
CA 842 42.70 43.07 -0.37 
CA 847 96.80 O26 -0.46 
CA 850 86.20 86.74 -0.54 
CA 849 64.80 65438 -0.58 
TABLE II 


In addition to the above comparisons, two runs were made between 
CH 258 and CH 259 as a further check on the latter to eliminate a 
possible source of error due to an incorrect value of gravity being 
used for the base station. The average of the two runs showed a 
difference of observed gravity between the two stations of 0.01 mgal 


less than the value published by Chapman [ 1966 ]}. 


We 
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Note: 
i) Plate I covers entire area 
ii) Plates IIa - Vf cover areas indicated 


Fig. 2 Diagram of Survey Area Showing 
Area Coverage by Photographic Plates 
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Fig. 11. Plate IIIb 
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V1i CONCLUSIONS 


—_———— 


The strong low SBA on Fort Ord indicates a depressed area with an 
isostatic root which 15 overcompensated In accordance with isostatic 
theory, it 1S expected that in the future this area should rise. As 
it rises the alluvium which presently covers the area will be removed. 
This process will continue until eventually the area is isostatically 
compensated. 

A detailed cross-section of the substructure was not attempted be- 
cause, as shown in Tables I and I], a reasonably accurate tie to the 


preceeding gravity work in the area could not be made. 


VII1. RECOMMENDATLONS 


“It is Strongly recommended that the entire Monterey Peninsula be 
resurveyed with a LaCoste-Romberg Mode! G Geodetic Gravity Meter or 
equivalent and the gravity anomalies be redeterminied accurately. 

Once the anomalies are well known a more thorough understanding of 
the subsurface structure can be obtained 

The only area on the west coast of the United States that has both 
the geology and negative anomalies consistent with the theory of sea- 
floor spreading is located basical!ly along the Salinas Valley. Thus 
Monterey is ideally located to study this concept. It is recommended 
that the areas both to the east of Fort Ord on ‘and and to the west 
at sea be surveyed to determine the full extent of the area of negative 
anomalies. This data in conjunction with a complete, accurate geologi- 
cal study of the same area should enable a full scale study of sea- 
floor spreading. 
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APPENDIX A 
Derivation of the Free Air Correction 
The normal gravity at any height h is expressed as: 
ch = 9, +{33 at (a) 


where gg, 9, equal normal gravity on the reference sphere and at 
elevation h, respectively. 


which expanded in a Taylor series gives: 


3 AY 2 5g, 3 
O,= 9 - (Shiloh + eGR), ne (1/6) 24) ts (b) 
which can be written: 
9, = Jo + Ah + Bht + Ch + Dh? +... (c) 


where A, B, C, etc. correspond to the coefficients of the Taylor 
expansion. 


Hirvonen [1960] developed the potential differences u = u, - U., 


in terms of a Taylor expansion along the normal and since 


ou. 

am 
and 

> ac 

0 hie ah 


his results can be applied here. 


Then 
2 2 Y¢ ee 2 ao 2 
A= kM [2+ e> - m- 27 me + e” - sin'L(3 & -5m-7— me 
7a 14 
a b 
+2e7 J 
and 
B= kM [3- 5m+ 2e° - sin? L, (9e° - 15m) ] 
2 2 
ab 


Ze 
where m= 4) a Db 
kM 
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k = universal gravitational constant 
M= mass of earth 

a = equitorial radius of ellipsoid 

b = polar radius of ellipsoid 


L = latitude 


e = second eccentricity of internationa! ellipsoid. 
When constants for the International ellipsoid are used 
a = 637838800 cm 
e*= 0.006722670022 


_m = 0,00344986385 cm3/gm - sec”? 





kM = 1.541363543 x 10°° sec ~* 

— 

A = - (0.308777237 - 0.000452063 sin® L, ) 

B = 07264891075 x 10° - 0.00021230891 x 10° sin? L 


ese 
Approximate the earth as a sphere of mass M and radius R, then 


considering only that part of the gravity due to the Newtonian at- 


traction 


the vertical rate of change of g with respect to R is 


dg _-2 kM 


which corresponds to the main term of A. 


Differentiation of the above yields 


d&g _6 kM 
7 4 
aR> er 


which corresponds to the main term of B. Also 


dg _ = 24 kM 
dR° R° 
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At any point Py on the earth's surface 


do KM 


Ro“ 


and thus 


3 
(Sa) 





g, = 980 gals 


Re =e Oo7 Ie king 
d3 

(2 | - -9,094392 x 10 2 gals/km° 
dR® O 


and 
2 -8 3 
te- - 1,532 x l0gs galsfkm se, 
By the same procedure 


0 4 
—z gals/km 


4 - 
(“4) aee0 20 y=) Gt 71371096 ee 10 
O 
dR Ro 
and 


-10 
D = 0.029737957 x 10. gals/km". 


Collecting terms and converting height to feet and gravity values to 
mgal: 


ye 
de> 95 (.0941154901 - .000137789 sin LOE 
-6 -6 .2 
= (.00678933167 x 10 - @°000T972422 xs10F Sim ESE 


Se 57oke0. 2 ee. (d) 


Truncating at 19° '° 


, the Free Air Correction is obtained from (d): 
FAC = (0.09411549 - 0.000137789 sine L)E 


- 0.0000000067 E¢ 
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